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The transverse momentum spectrum of rj meson in relativistic heavy-ion collisions is studied at 
the Next-to-Leading Order (NLO) within the perturbative QCD, where the jet quenching effect 
in the QGP is incorporated with the effectively medium-modified rj fragmentation functions using 
the higher-twist approach. We show that the theoretical simulations could give nice descriptions of 
PHENIX data on rj meson in both p-|-p and central Au-|-Au collisions at the RHIC, and also provide 
numerical predictions of rj spectra in central Pb-|-Pb collisions with ^sjvat = 2.76 TeV at the LHC. 
The ratios of 77/71° in p -|- p and in central An -|- An collisions at 200 GeV are found to overlap in 
a wide pr region, which matches well the measured ratio 77/71° by PHENIX. We demonstrate that, 
at the asymptotic region when pr —>■ 00 the ratios of 77/77° in both An -|- An and p -b p are almost 
determined only by quark jets fragmentation and thus approach to the one in e''"e“ scattering; in 
addition, the almost identical gluon (quark) contribution fractions to 77 and to tt result in a rather 
moderate variation of 77/77° distribution at intermediate and high pr region in A-|-A relative to that 
in p+p; while a slightly higher 77/77° at small pr in Au-|-Au can be observed due to larger suppression 
of gluon contribution fraction to 77° as compared to the one to 77. The theoretical prediction for 
77/77° at the LHC has also been presented. 
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The strong suppression of single hadron production at 
large transverse momentum [l|, Q has provided the con¬ 
vincing evidence of the jet quenching phenomena discov¬ 
ered in relativistic heavy-ion collisions (HIC) Q. Ex¬ 
tensive phenomenological investigations and experi¬ 
mental measurements dMli on the suppression of single 
hadron spectra at high pp have been carried out at both 
the RHIC and the LHC. As the first observable of jet 
quenching phenomena, the yield suppression of inclusive 
hadrons is arguably the most thoroughly studied quan¬ 
tity of jet quenching, and provides an indispensible tool 
to extract the properties of the hot medium created in 
nucleus-nucleus collisions by comparing theoretical calcu¬ 
lations with experimental measurements, such as the jet 
transport coefficient q 0. The interplay between the¬ 
ory and experiment on the single hadron production will 
help constraining the longitudinal distribution of parton 
energy loss in hot/dense QCD medium, and better un¬ 
derstanding the jet-medium interactions after being com¬ 
bined with studies of full jets which also shed light on the 
angular distribution of the medium-induced gluon radi¬ 
ation and thus constrain the transverse distribution of 
parton energy loss as well 18-^. 


So far, most of the theoretical calculations on single 
hadron productions in HIC focus on tt meson or charged 
hadrons (where tt also giving a predominant contribu¬ 
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tion), and there are very few studies on other identified 
hadrons 24 -2^. We note that 77 meson is the second 
important source of decay electrons and photons just af¬ 
ter the 77°, and its quantitative analysis is of importance 
to suppress the noise in measurements such as direct 
photon [2^ [^. In addition, PHENIX has made a de¬ 
tailed measurement of 77 production in Au-|-Au collision 
at ^/sNN = 200 GeV, but there is not any theoretical 
calculations on 77 at high px in HIC to the best of our 
knowledge. It is of great interest to see how parton en¬ 
ergy loss effect alters 77 spectrum in HIC and whether the 
theoretical model of jet quenching could make a simulta¬ 
neous description of both tt and 77 production. 


In this Letter we study 77 productions at large px in 
high-energy nuclear collisions at the RHIC and LHC in 
the framework of higher twist approach of jet quench¬ 
ing In higher twist approach, multiple scatter¬ 

ing of the fast parton traversing through the dense QCD 
matter is calculated with generalized QCD factorization 
of twist-4 processes. The resulting parton energy loss 
due to multiple scattering leads to effectively modified 
parton fragmentation functions in medium (mFF). In¬ 
corporating mFFs into a next-to-leading order (NLO) 
perturbative QCD (pQCD) improved parton model, a 
phenomenological study on 77° suppression in nucleus- 
nucleus collisions [1, gives a decent description of 77° 
yields in Au-|-Au collisions at the RHIC and in Pb-|-Pb 
reactions at the LHC. In this Letter, we may employ the 
identical model to investigate px spectrum of 77 produc¬ 
tion at NLO in HIC, with the same jet transport param- 
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FIG. 1 : Comparison between the PHENIX data of rj invari¬ 
ance cross section in p + p collisions at 200 GeV and the NLO 
pQCD theoretical calculations. 



FIG. 2: Comparison between the PHENIX Data of rj nuclear 
modification factor in Au -I- An collisions at 200 GeV and nu¬ 
merical simulations at NLO. 


eters extracted in 7r° production in HIC. We also explore 
the features of ij/rT^ ratios in both p -b p and A -b A colli¬ 
sions. Although the flavor dependent parton energy loss 
in the QGP may alter the flavor compositions of fast par- 
tons which in turn have distinct probabilities fragment¬ 
ing to rj and to tt, at very high px, the predomination 
of quark contribution for both rj and 7r° ensures that the 
QCD medium effect does not affect the productions ratio 
of rj and tt®. The almost identical gluon (quark) contri¬ 
bution fractions to rj and to tt in p -b p may lead to a 
rather small deviation of rj/r:^ in A -b A with respect to 
that in p -b p at intermediate and high px region. And 
at small px, an enhancement of rj/n^ in A -b A collisions 
can be seen because of the relative larger suppression of 
gluon contribution fraction to tt® than that to 77 in A -b A 
collisions. 


Firstly, we investigate the single inclusive rj produc¬ 
tion in p -b p collisions, which provides the baseline for 
the nuclear modification of 77 yield in HIC. A NLO pQCD 
improved parton model for the initial jet production spec¬ 
tra has been employed 3^ , in which the inclusive particle 
production cross section in p -b p collisions can be factor¬ 
ized into a convolution of parton distribution functions 
(PDFs) inside the proton, elementary partonic scattering 
cross sections da/dt, and parton fragmentation functions 
(FFs), 


da: 


,72 — ^ ^ / ^^adxi}fa/p{Xa-,(^ ) 

^ ^ abed ^ 


x—j{ab^cd) - \-0{ag), ( 1 ) 

di 


where dd-/dt{ab —>• cd) denotes leading order (LO) par- 


ton scattering cross sections at a^. The NLO partonic 
cross section at at includes 2 —>■ 3 tree level contri¬ 
butions and one loop virtual corrections to 2 —> 2 tree 
processes [S^. In the above equation fa/p{xa, stands 
for the PDFs in proton with Xa the momentum fraction 
of the beam proton carried by the incoming parton, and 
in numerical simulations we use CTEQ6M parametriza- 
tion [3^. D^^^{zc, represents the parton fragmen¬ 

tation function in vacuum, which gives the possibility 
of parton c fragmenting into hadron h with momentum 
fraction Zc- The factorization scale, renormalization scale 
and fragmentation scale are usually chosen to be the same 
and are related to px of the final hadron. This pQCD im¬ 
proved parton model at NLO accuracy has worked very 
nicely to describe data on inclusive pion production in 
p-bp collisions at RHIC by using AKK pion FFs and the 
scale in the range /r = 0.5 ~ 1 . 5 px [^.13^. 

To compute rj production at LO and NLO with pQCD, 
the non-perturbative input of 7 ] FFs will be needed. The 
availability of a parametrization of rj FFs AESSS [s^ 
allows us to make a NLO pQCD calculations for single- 
inclusive eta meson production as a function of final state 
px in hadron-hadron collisions. In AESSS parametriza¬ 
tion because of the absence of enough data on inclusive 
77 productions, the 77 FFs can not be extracted separately 
for each quark flavor without additional assumptions, 
and the assumption is made that all light quark frag¬ 
mentations are the same. In Fig . [J we confront our cal¬ 
culation with PHENIX data [39j, and it is observed that 
the computed inclusive 77 spectrum in p+p collisions with 
the scale ji = l.Qpx agrees well with the PHENIX data. 


In central high-energy nucleus-nucleus collisions, a hot 
and dense QCD matter may be created. A fast parton 
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FIG. 3: Predictions of rj medium modification factor in FIG. 4: Comparison between the PHENIX data of > 7 / 71 ° ratio 
Pb + Pb collisions at 2.76 TeV with different qo parameters. inp + p collisions and An + An collisions at 200 GeV and the 

numerical simulations at NLO. 


propagating in the QCD medium may encounter multiple 
scattering with other partons in medium and lose energy 
by induced gluon emission. In higher twist approach of 
jet quenching this kind of multiple scattering is 

related to twist -4 processes of hard scattering and may 
give rise to the effective modifications of parton FFs in 
vacuum such as: 


D'^{zh,Q^) = + 


2 tt 


3 ^ 


dz 




+ l^lq^gq[z,X,XLjT)D'l[ — ,Q‘^) ,( 2 ) 


where A^q^qg{z,x,XLA\) and Ajq^gq{z,x,XL,iT) = 
A7g_>gg(l — z,x,xl,£‘^) are the medium modified split¬ 
ting functions calculated in higher twist approach [s^, 
[^ . The medium-modified FFs averaged over the ini¬ 
tial production position and jet propagation direction are 
given as [1,0: 


{D'^izh,Q^E,b)) = 


J d^rtAiM)tBi\b - T^) 


d(j) 

27 r 


d^rtA{\r\)tB{\b - r\)b'^{zh, E, r, (j), b). 

( 3 ) 


where b is the impact parameter and Ia b are the nuclear 
thickness functions in Glauber model jdl]. In the higher 
twist approach, one assumes that the parent energetic 
parton loses its energy in the QCD medium prior to vac¬ 
uum hadronizations and single hadrons at large px are 
fragmented in vacuum by the partons passed through the 
QCD medium. 


Therefore, cross section of the single hadron in HIC 
collisions could be expressed as: 


1 


da'^ 


AB 


jAB(ia ^ 

abed '' 


^hinib) dyd?pT 


dXadXbfa/AiXa,P^)fb/B{xb,fb 


X -^(ab -A cd) 

dV ’ 


{b’^izh,Q^E,b)) 


nze 


+ 0{al). 


( 4 ) 


Here (b) gives the number of binary nucleon-nucleon 
collisions at the impact parameter b in A + B collisions, 
fa/A{xa, pb represents the effective PDFs inside a nu¬ 
cleus. In our calculations, we employed EPS 09 NLO nu¬ 
clear PDFs to include initial-state cold nuclear matter 
effects 411 . 


We assume all the energy loss of a fast parton is that 
carried away by the radiative gluon in the multiple scat¬ 
tering processes, the corresponding parton energy loss in 
the QCD medium can be expressed as: 


AE 


N, 


dy dzdi^ 


(1 + Z )3 


'xqR{E,y) sin^ 


j 


T 

T 


AEx{l - z)‘ 


( 5 ) 


which is proportional to the jet transport parameter 
qB.iE,y). The jet transport parameter is related to the 
parton density distribution in the medium, therefore the 
space-time profile of the jet transport parameter char¬ 
acterizes the medium properties in the medium modi¬ 
fied fragmentation functions and energy loss calculation. 
Phenomenologically, we treat it proportional to the par- 
ton density relative to their values at the center of over¬ 
lapped region in the most central collisions at initial time 
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FIG. 5: The NLO pQCD theoretical prediction of 77 /vr^ ratio 
in Pb + Pb collisions at 2.76 TeV. 


and a given temperature in an ideal hydro or in hadronic 
phase, the four momentum of the jet and the four flow 
velocity in the collision frame along the jet propagation 
path are also included 

A full three-dimensional(3+lD)ideal hydrodynam¬ 
ics [i^, 431 to describe the space-time evolution of the 


QCD medium in heavy-ion collisions is employed to give 
the required space-time evolutionary informations of the 
medium such as: parton density, temperature, fraction of 
the hadronic phase, four flow velocity. In this theoretical 
framework, there left only one parameter when describe 
the medium modification of the hadron productions in 
heavy ion collisions: qoTQ. In this Letter, we fix the ini¬ 
tial time of the QGP medium at tq = 0.6 fm, and de¬ 
termine the initial values of the jet transport parameter 
qo by fitting the final charged hadron multiplicity den¬ 
sity in mid-rapidity of the hydrodynamics results [ 1 , . 

The qo can also be considered to adjust the strength of 
jet-medium interaction as a parameter, the larger qo is, 
the stronger the jet-medium interaction will be at every 
space-time points. 


With all these ingredients, we can calculate the hadron 
productions in heavy ion collisions up to the NLO. Al¬ 
though we can directly compare the calculated hadron 
productions in heavy ion collisions with the experimen¬ 
tal data, we introduce the nuclear modification factor 
Raa to better demonstrate and understand the medium 
modification phenomenon: 


n ^ da’XB/dycfpT 


( 6 ) 


At various values of qoTo = 0.48 — 1.02 GeV^,we calcu¬ 
late the Au + Au productions within medium modified 77 


fragmentation functions in the 5% most central Au + Au 
collisions at RHIC energy -/s = 200 GeV, and compared 
our results with the PHENIX experimental data on Raa 
of the 7] spectra. The theoretical calculation results ex¬ 
plained the data well in large pr region in Fig. [21 The well 
agreement with the PHENIX data 28| at mid-rapidity in 
the range pt = 2 — 20 GeV shows that, even 77 meson is 
4 times heavier than 7 r°, a similar flat production sup¬ 
pression has been observed at RHIG in this px range 
independent of their mass. 

The prediction of Raa of the 77 spectra in the central 
Pb + Pb collisions at the LHC energy = 2.76 TeV is 
given in Fig. [3] The go'^o are chosen from 0.84—1.8 GeV^. 
The suppression factor increases with px was observed, it 
is mainly due to the energy dependence of parton en ergy 
loss and the less steep initial jet production spectra [4^ . 

The px dependence of the ’q/iT^ ratios are also investi¬ 
gated. We plot the px dependence of the 77 / 71 ° ratios in 
Au + Au at 200 GeV as in Fig. |4l and a good agreement 
between the model calculations with PHENIX data can 
be seen. We also predict the px dependence of the 77 / 71 ° 
ratios in Pb -I- Pb at 2.76 TeV in Fig. |5l Similar trend 
could be seen at the RHIC and LHC that with the in¬ 
creasing of Px, the 77 / 77 ° ratio in A -|- A collisions comes 
closer to the p+p curve, and at very larger px two curves 
coincide with each other. 

Due to the different color factors of quark-quark vertex 
and gluon-gluon vertex, in QCD medium gluon jet suf¬ 
fers larger energy loss than quark jet. Because quark FFs 
and gluon FF into 77 and tt have quite different features 
(for instance, the bottom panel of Fig. [ 6 ] shows that at 
very highpT region, Dg_>,)/Dg_>,^o at z/i = 0.7 is approx¬ 
imately 0.5), in principle, a change of flavor compositions 
of parton jets may affect the ratio of 77 / 77 ° [ 2 ^. The coin¬ 
cidence of the A -I- A and the p -I- p 77 / 77 ° curves in a wide 
Px region can not be explained in one simple story that 
parton jets loss their energies hrst in the QCD medium 
and then fragment into hadrons in the vacuum [28l |. 


To have a better understanding of the ratio of differ¬ 
ent hadron yields, we simplify the formula of identified 
hadron px yield in p -|- p collisions as: 


1 dcr^o 
Pt dpx 


= j fqi^) ■ Dq^ri.7vo{Zh,Px) 

[ fg{ — )-Dg^ri,^o{zh,Px)- 
J Zh ^ ” Z 


dzh 


(7) 


The above equation shows that the hadron yield in p -|- p 
will be determined by two factors: the initial (parton-)jet 
spectrum fq^g{px) and the parton fragmentation func¬ 
tions Dqg^gT^o{zh,px)- We plot the gluon’s portion in 
total parton FFs to 77 and to 77 ° at fixed z = 0.7 in 
Fig. IlKtop panel), we can see that, at very high px, 
(up-)quark FF is much larger than gluon FF for both 77 
and 77° and shows weak px dependence. At very high 
Px, the yields of 77 ° and 77 will be almost only given by 
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FIG. 6: Top panel: gluon portion in total FFs (top panel) at 
fixed z = 0.7; Bottom panel: ratio of the gluon, quark 

fragmentation functions (bottom panel) at fixed z = 0.7. 


quark fragmentation (also see top panel of Fig. [7|), and 
we can neglect the contributions of gluons. Noticing that, 
at high pt, quark FFs Dq^^ Tro(zh, Q = pr) have a weak 
dependence on Zh and pr in the typical Zh region 0.4 —0.7 
for identified hadron production |45| . In the asymptotic 
region with oo we obtain 


R{p/n°) 


dar, / 
dpT 


ddT^O 

dpT 


j fq{^).Dq^^o{Zn,PT)^ 
SgFlq—,P t) 
^qDq^T,o{{zh) ,Pt) 


( 8 ) 


Therefore, even though quarks and gluons may lose 
different amount of energies and the flavor content of 
parton-jets will be changed in A + A collisions, at the 
very high px region in A + A collisions, the yields of both 
TT^ and p should also predominantly come from quarks, 
and the ratio p/ir^ will also be determined only by quark 
FFs in vacuum with Zh shift due to energy loss effect. 
As we mentioned before, in this kinematic region quark 
FFs have a weak dependence on Zh and px, so we reach 
the conclusion that, at very high px region, the ratios of 
p/tt^ in both A +A and p + p should overlap with the one 
in e“*"e“ scattering, and reach a universal value ^ 0.5. 


Next, we consider the feature of p/ir'^ when the trans¬ 
verse momentum px is not very high. We define 


(px) = 


1 


Pt dpT 


(9) 


which denotes the gluon contribution fraction to p and 
TT^ yields with Eq. 0. We plot the gluon (quark) con¬ 
tribution fraction to p and yields in p -|- p collision in 
Fig. [7] (top panel). It shows that the quark contribution 
indeed dominate the p and 7r° productions at larger px 
region. Surprisingly, we also observe that 

G^\px)-G^{px) 


in the px region of 4 — 20 GeV in p -|- p at RHIC. 
From Eq. 0 we derive 


R{p/tt°) = 


Jfq{P^).Dq^,{z,,PT)^ 
Ifqi^)-Dq^MZH,PT)^ 
Ifqi^^)-Dq^,{zu,Px)^ 
^f9m-Dq^.o{zn,PT)^ ■ 


( 10 ) 


It implies that whereas and p yields in p-|-p come from 
quark or gluon hadronization, when calculating the ratio 
ry/TT^, because of the relative identical fractional contri¬ 
butions of gluon and quark to ir^ and ry, we can consider 
the contributions of only quarks (or gluons). The flavor 
compositions or mixture of quarks and gluons in p -|- p 
have nearly negligible effect on p/tt®, which is confirmed 
by the numerical results shown in the bottom panel of 

Fig.m 


Now we turn to p/tt^ in A+A collisions. Though flavor 
dependent energy loss in hot QCD medium may change 
the flavor compositions of jets, i.e., gluons may suffer 
more suppression than quarks, the effect to p/tt^ will be 
minimized by the nearly identical G'^ {pt) and G'^{px) 
observed in p-l-p, which underlies the overlapping of p/tt^ 
in A -|- A and the one in p -|- p in a wide region of px. 

Nevertheless, small modifications to p/7r° due to jet 
quenching effect will still exist. In A -|- A the relation 
between (pr) and G'^{px) will be broken a little bit. 
We plot the gluon and quark contribution fractions to p 
and 7r° yields in Au -|- Au collisions at 200 GeV in Fig. [S] 
It reveals that the flavor dependent jet quenching effect 
really results in smaller gluon contribution fraction in 
heavy-ion collisions as compared to p -I- p in Fig. 0 A 
naive expectation is that because gluon may give larger 
p/tt® ratio than quark does (see bottom panel of Fig. [6|), 
the larger suppression of gluons in the QCD medium will 
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FIG. 7: Top panel: quark and gluon contribution fractions 
to total 7 r° (or rj meson) yields at NLO in p + p collisions at 
200 GeV. Bottom panel: the ratio rj/n^ at NLO when only 
gluon or up quark contribution is considered in p + p collisions 
at 200 GeV. 


reduce rj/ir^. However, contrary to this naive expectation 
the jet quenching effect may instead enhance the ratio 
rj/ir^ in high-energy nucleus-nucleus collisions as seen in 
Fig. m and Fig. [5] The reason is that the suppression 
of gluon in QCD medium imposes a larger reduction of 
the yield of 7r° than that of rj (top panel of Fig. |S]), thus 
giving rise to a slightly enhanced rj/ir'^ ratio in A-|-A. To 
demonstrate this we plot in the bottom panel of Fig. [8] 
the rj/n^ ratio in heavy-ion collisions originated from only 
gluon or quark contribution. One can see that a larger 
rj/ir^ ratio with only gluons included in the theoretical 
simulations. We emphasize that the identified hadron 
yield in heavy-ion collisions relies on three factors: the 
initial hard jet spectrum, the energy loss mechanism, and 
parton fragmentation functions to the hadron in vacuum. 
When rj/ir^ ratios or other ratios of hadron yields are 
concerned , a simple estimation with only last two factors 
such as (flavor dependent) energy loss effect and FFs in 
vacuum may sometimes reach a wrong conclusion. 
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FIG. 8 : Top panel: quark and gluon contribution fractions to 
total 7 r° (or 77 meson ) yields at NLO in Au -I- An collisions 
at 200 GeV. Bottom panel: the ratio rj/'K^ at NLO when 
only gluon or up quark contribution is considered in Au -|- Au 
collisions at 200 GeV. 
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